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Abstract 

Simulations that couple different classical molecular models in an adaptive way by changing 
the number of degrees of freedom on the fly, are available within reasonably consistent theoretical 
frameworks. The same does not occur when it comes to classical-quantum adaptivity. The main 
reason for this is the difficulty in describing a continuous transition between the two different kind 
of physical principles: probabilistic for the quantum and deterministic for the classical. Here we 
report the basic principles of an algorithm that allows for a continuous and smooth transition by 
employing the path integral description of atoms. 
PACS numbers: 02.70.Ns, 61.20.Ja, 61.25.Em 



Introduction: The development of adaptive resolution simulation schemes is a subject of 
growing interest within the community of condensed matter, material and chemical physics. 
By adaptive resolution it is meant that the space is partitioned in regions characterized by 
different molecular resolutions where molecules can freely diffuse changing their representa- 
tion according to the region where they are instantaneously located. In the last few years 
several approaches have been presented and they are characterized by different levels of 
theoretical sophistication and computational complexity jl-6]. The interest in this kind of 
approach arises from the fact that it may efficiently tackle the problem of interplay between 
different scales. However while the adaptive process can be described in a reasonable way 
according to the basic principles of classical dynamics and thermodynamics, the same can- 
not be said when quantum mechanics enters into the game. In general, the proper coupling 
of quantum and classical mechanics is known to be a non trivial (and open) problem (sec 
e.g.(7|) and here the adaptive character adds up as a further major difficulty [8j. In this work 
we aim to develop an approach where the idea of coupling a classical and quantum molecular 
model in an adaptive fashion can be made in a way that the "probabilistic-deterministic" 
discontinuity is removed and yet the adaptivity takes place in a smooth controlled numerical 
way. 

Path Integral representation of atoms and molecules: In this work our quantum sys- 
tems are composed by atoms described within the path integral formalism. This approach 
is, by now, a standard tool in molecular dynamics and well established in literature (see e.g. 
[9|). For this reason here we do not present a technical description of the method, instead we 
explain why this approach is optimal to adaptively couple quantum systems with a classical 
bath. The important aspect of the path integral approach for this work is that an atom, 
which is usually represented as a sphere in classical force fields, becomes, within the path 
integral formalism, a classical "polymer ring" so that the interaction site is delocalized into 
the beads of such a polymer; each bead of the polymer is linked to its next neighbors along 
the chain by a harmonic potential. The elastic constant depends on the temperature, T, 
of the system and the number of beads, n, used to represent each atom, k = m ^?iP ; with 
&b being , as usual, the Boltzmann constant, m the mass of the atom, h Planck's constant. 
Such a derealization characterizes the quantum nature of the atom. This approach has been 
extensively used in molecular dynamics to account for basic quantum effects otherwise not 
described by standard classical force fields 10 is| . For the adaptive process the classical 
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polymer representation of atoms has far reaching consequences, because it translates the 
quantum-classical coupling into the coupling of two effective classical regions characterized 
by a different number of (as a matter of fact) "classical" degrees of freedom; thus the whole 
machinery of classical adaptive methods would apply straightforwardly. In this work we 
show that indeed this is the case. Obviously the approach here is suited only for a class of 
problems and, as it is designed now, certainly not for those problems where the quantum me- 
chanics refers to electrons. It must be noted that practical methods H] which couple 
the two regimes when electrons are considered suffer from the same conceptual limitations 
underlined in this work and base their validity on numerical criteria only. In this specific 
case the conceptual discontinuity is in the arbitrary cut off of the electron wavefunction and 
in the non conservation of the number of electrons as the system evolves in time. Although 
we do not consider systems with electrons, the concept of mapping a quantum problem onto 
a classical one and then have a smooth adaptivity is anyway of general valence. In this sense 
this work represents a possible starting point for a general adaptive coupling of quantum 
and classical description of atoms and molecules. 

The Adaptive scheme: Regarding the adaptive scheme for the classical case, the basic 
requirement is a controlled procedure of changing the number of degrees of freedom based on 
solid physical principles and consistent with the thermodynamic equilibrium of the overall 
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4201] . In this respect the AdResS method meets the requirements above in a (ex- 



tensively tested) satisfactory way 2jJ-|25j. For this reason in this work the coupling between 
the polymer rings and the classical particles will be done within the AdResS framework. Be- 
low we introduce the basic features of AdResS relevant for this work. In the classical AdResS 
method the atomistic (high resolution) and the coarse-grained (low resolution) regime are 



coupled via a force interpolation (see e.g. 

F a/3 = w(X a )w(X fi )Fy + [1 - w(X a )w(X p )]F^ (1) 

where a and indicates two molecules, F atom is the force derived from the atomistic force 
field and F cg from the corresponding coarse-grained potential, X is the x coordinate of the 
center of mass of the molecule and w is an interpolating function which smoothly goes from 
to 1 (or vice versa) in a transition region (A) where the lower resolution is then slowly 
transformed (according to w) in the high resolution (or vice versa), as illustrated in FigJT] 



(see also 



26j). An additional locally acting thermostat is employed to assure the overall 



thermodynamic equilibrium. The coarse-grained potential is obtained from a reference all 
atom simulation at the given thermodynamic condition via an iterative inverse Boltzmann 
procedure employing the molecular center of mass radial distribution [27J]. For the case of 
the path integral description we have exactly the same procedure as for the classical case 
but with the difference that instead of a classical atomistic representation of the molecule 
we have a molecular representation where the atoms are described as polymer rings in a 
path integral approach (see FigJT]). It follows that the coarse-grained model is derived from 
a full path integral reference simulation at the given thermodynamic condition. Given the 
framework reported above the coupling between the path integral representation and the 
coarse-grained model occurs via the interpolation according to w between F pi acting on the 
beads of the rings (which now plays the equivalent role of Y atom in EqffJ) and the F° 9 derived 
from the coarse grained potential acting among the centers of mass of the molecules. We 
tested this idea studying a liquid of tetrahedral molecules whose atomistic model was used 
in the original development of AdResS (see FigJU(b)). 

Results of the Path integral-coarse-grained Adaptive Resolution Simulation: We 

studied a system of thousand molecules at two different temperatures, indicated as T\ and 
T 2 . T\ is the same temperature employed in the previously studied classical model and 
T 2 = -^=. In the Lennard- Jones units of this paper (see the Appendix or Ref. |21|) T\ — 1 



and consequently T 2 = -7=. The different temperatures are directly related to the elastic 



constant of the polymers and thus, for the same number of beads per ring (n = 10), a 
lower temperature means more flexibility of the polymer rings, thus we have k% = For 
the testing purpose we have deliberately chosen T 2 because it mimics the thermodynamic 
conditions of a "more quantum" system than T\ (see also |28|). Given the computational 
cost of path integral simulations and the very extended range of tests we have used, for 
both systems, the number of beads per polymer, n, equal to 10 29[]. Figj2] reports the radial 
distribution functions (RDF) and the density distribution of the AdResS simulation of thou- 
sand molecules compared to a full path integral simulation for the temperature T\. They 
show that indeed this coupling procedure displays the desired behaviour. In particular, in 
FigfJ^b), the comparison between the bead-bead radial distribution function obtained with 
AdResS in the quantum region and that obtained from a full path integral simulation shows 
that indeed the very quantum nature of the particles in the quantum region of AdResS is 
very well described and thus the "classical bath" of the coarse-grained molecules is indeed 
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able to reproduce the overall thermodynamic conditions properly. Regarding the density 
in Figj2]^c), the deviations with respect to the reference value (above all in the transition 
region) are the same obtained in the classical case 21] . As for the classical AdResS case, 
one must be sure that indeed there is exchange of molecules between the two regimes. In 
fact it may occur situations where the two regimes are in equilibrium because some barriers, 
artifact of the method, would hinder the free exchange so that the molecules are reflected 
back at the transition region. Figj3] shows that this is not the case and the molecules diffuse 
between the two regions in a proper way. FigJH shows the results for T 2 , and also in this case 
the agreement with the reference full path integral simulation is satisfactory (the diffusion 
profile, not reported here, shows also the same behaviour as for the case Ti). The overall 
agreement make us confident that the principles of the adaptive algorithm can be extended 
to the case of path integral representation. 

Conclusions: In general for physical systems one should not expect the agreement found 
here although it can be expected a reasonable accuracy in the same way (and with similar 
limitations) of the classical AdResS. Having shown that the idea works in a rather reasonable 
way, it remains to address the question of why such an approach may be useful for appli- 
cations and not a mere conceptual exercise. In general, there are several problems in soft 
condensed matter where relevant properties are the results of the interconnection between 
the local bonding and the larger scale molecular packing. In such cases classical models may 
give a satisfactory description of the large scale packing but are not sufficient for the high 
resolution required locally and the path integral description becomes crucial ; this is the 
situation where our approach would represent an optimal tool. Another important field of 
application is that involving liquid water and above all water as a solvent of large molecules 
\3l\ . In particular in cases where the solute- water hydrogen bond plays a crucial role 32|-|34| 
and classical model can only partially describe the water mediated process. Moreover, com- 
bined with specific force fields [35], would be ideal to study local protonation/deprotonation 
of molecules in solution. In these cases the adaptive method would require that only the 
region around the solute should be treated with the path integral resolution while the rest 
can be treated in a coarse-grained fashion as shown here or with a hierarchy of models as 



already done in part within AdResS 2% |25|, |36] . Needless to say that this is only a first step 



towards the development of a generic scheme for quantum-classical coupling and serve as a 
basis for further work. 
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Appendix: The Classical atomistic model 

Here we briefly report the basic technical details of the force field for the tetrahedral 
molecules. Each molecule is characterized by N = 4 atoms of the same mass m connected 
by anharmonic bonds. All atoms in the system interact according to a purely repulsive 
shifted 12-6 Lennard- Jones potential with a cutoff at 2 1//6 cr: 

f Ae\(^-) 12 - (^H 6 + ±1; r iaj0 < 2 1 /6 a 
[ 0; r iaj p > 2 1/6 cr 

fiajp defines the distance between the atom ia of the molecule a and the atom j/3 of the 
molecule 0. We define e as a unit of energy. All atoms have the same excluded volume 
diameter a, where a is the unit of length. Neighboring atoms of a molecule are linked via 
an attractive FENE potential 

rratom , ,1 4^0 Ml ~ (^fll T iaja < R 

u FENEV iaja) — \ W) 

I oo; r iaja > R 

Rq = 1.5a is the divergence length and k = 30e/a 2 the stiffness. 
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FIG. 1: Pictorial representation of the adaptive box and molecular representation. The region 
on the left, indicated by CG, is the low resolution region (coarse-grained), the central part is the 
transition (hybrid) region A, where the switching function w{x) is defined, and the region on the 
right, indicated by PI is the high resolution region (full path integral region). 
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FIG. 2: Results for the adaptive simulation of the liquid of tetrahedral molecules at the temperature 
T\. Top (a), the molecular center of mass-center of mass radial distribution function obtained with 
AdResS is compared with that obtained from the full path integral reference system,. Middle (b), 
the bead-bead radial distribution function obtained with AdResS in the quantum region compared 
with that of the full path integral reference system. Bottom (c), the particle density in AdResS 
compared with the reference system. The density is equal to 0.1<r~ 3 in the units reported in the 
Appendix. 
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FIG. 3: Diffusion profile for the molecules moving from the path integral region to the coarse 
grained one and for molecules moving in the opposite direction, for the system at temperature T\. 
The picture shows that no barriers, due to possible artifacts of the algorithm, hinter the diffusion 
process. 
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